Takada et al. SpringerPlus 2014, 3:418 
http://www.springerplus.eom/content/3/1/418 



0 SpringerPlus 

a SpringerOpen Journal 



RESEARCH 



Open Access 



Comparison of generalist predators in winter- 
flooded and conventionally managed rice paddies 
and identification of their limiting factors 

Mayura B Takada^ , Shun Takagi^, Shigeki Iwabuchi^ Takuya Mineta"^ and Izumi Washitani^ 



Abstract 

Winter-flooding of rice paddies witliout tlie application of agricultural chemicals is attracting attention as a 
new agricultural method for enhancing the habitat conditions of wintering waterfowl in rice paddy ecosystems 
throughout Japan and east Asia. Conditions in these paddies are expected to result in restoration of not only the 
winter habitats of waterfowl but also those of other taxonomic groups during the rice growing season. In this 
study, we tested whether the diversity of summer spiders — ubiquitous predators in rice paddies — was higher in 
the winter-flooded paddies than in the conventional ones by conducting field measurements in 31 winter-flooded 
and 7 conventional paddies. Limiting factors of spiders in the winter-flooded paddies were then examined. Results 
revealed that both the density and species richness of spiders were significantly higher in the winter-flooded paddies than 
in the conventional ones both before and after the insecticide application against pecky rice bug Stenotus rubrovittatus 
(Matsumura)(Hemiptera: Miridae) to conventional paddies. In addition, spider density and species richness in the 
winter-flooded paddies correlated with the availability of two prey groups — chironomids and other nematocera. 
These findings suggest that in the winter-flooded paddies the diversity of generalist predators is higher than in 
the conventional ones during the rice-growing season and that the combination of management at both the 
landscape and field level is likely more effective for increasing spider abundance in winter-flooded paddies. 
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Introduction 

Biological communities in agricultural crop fields that 
have received low or no chemical application are more 
diverse than those of conventionally farmed fields (e.g. 
Bengtsson et al. 2005). In particular, rice paddies main- 
tained by biodiversity-friendly farming techniques are 
expected to function as important substitute habitats for 
organisms in natural wetlands, which have declined due 
to the development of urban and agricultural fields (Fasola 
and Ruiz 1997; Washitani 2007). 

The farming practices of biodiversity-friendly paddy 
fields include the winter flooding of paddies to enhance 
the habitat conditions for wintering waterfowl (Elphic 
and Oring 1998; Tourenq et al. 2003; Kurechi 2007). In 
the paddies at the Ramsar wetland site that surrounds 
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the natural Kabukuri-numa Wetlands, which is one of 
the most important wintering sites for geese in Japan, 
the winter-flooding regime combined with no agricul- 
tural chemical applications has become a growing trend. 
This farming method is attracting attention as a new 
biodiversity-friendly approach throughout Japan and 
East Asia (Kurechi 2007; Washitani 2007). In fact, to re- 
store paddy biodiversity including wintering waterfowl, 
no agricultural chemicals have been applied to the winter- 
flooded (WF) paddies in this region since around 1998. 
Moreover, WF paddies are also expected to enhance eco- 
system services for rice production, including the biological 
control of insect pests by generalist predators. 

In this study, the diversity of spiders — ubiquitous 
predators in rice paddies — was compared between WF 
and conventional paddies through field measurements 
conducted in 31 WF paddies and 7 conventional ones. 
In addition, factors limiting the abundance of spider 
populations in WF paddies were also examined. Crop 
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yield in biodiversity-friendly fields can be low com- 
pared with that in paddies that receive conventional 
chemical applications (de Ponti et al. 2012). Thus, 
effective pest management techniques other than the 
use of insecticides are needed to sustain growth in 
biodiversity-friendly farming (Zehnder et al. 2007). Identify- 
ing the limiting factors of generalist predators in the WF 
paddies can help with the construction of effective 
management plans for insect pest control, because 
our previous studies showed that the spiders predate 
insect pests (Kobayashi et al. 2011; Takada et al. 
2013), and suppress the density and the consequent 
damage in organic paddy fields (Takada et al. 2012). 
Spider population growth is known to be constrained 
by habitat complexity (Langellotto and Denno 2004; 
Takada and Miyashita 2004, 2014) and/or prey 
availability (Wise 1979; Miyashita 1992). In this 
study, we focused on the vegetation structure of rice 
and abundant weeds, as well as the biomass of abun- 
dant arthropod groups, such as chironomidae and 
other nematocera, brachycera, and hoppers including 
Cicadellidae and Delphacidae, and the pecky rice bug 
Stenotus rubrovittatus, as potential limiting factors of 
the spiders. 

Methods 

Study site 

A field survey was conducted in the area surrounding 
Kabukuri-numa at Osaki City, Miyagi Prefecture, Northern 
Japan (38°37'N, 141°07'E), where many WF paddies (total 
area: about 0.2 km^) have been cultivated (Kurechi 2007) 
since around 1998. The mean annual precipitation from 
2002 to 2007 was approximately 1166 mm, and the mean 
temperature in August during the same period was 23.6°C 
(Japan Meteorological Agency 2012). In this region, 
insecticides for Stenotus rubrovittatus (Hemiptera: 
Miridae), the most abundant and influential rice pest 
(Yoshioka et al. 2011; Takada et al. 2012; Yoshioka 
et al. in press), are applied to most of the conventional 
paddy fields after heading of rice plants in early 
August, when spider diversity in paddies is at the high- 
est level (Oyama et al. 2005). Field measurements 
were conducted twice, just before and after application 
of the insecticides to conventional fields, with the 
expectation that spider diversity decreases in conven- 
tional paddies, while remaining unchanged in WF ones 
after the application of insecticide to conventional 
fields. 

In this area, clusters of WF paddies were distributed 
among conventionally managed paddies. The WF fields 
have been maintained by using similar management 
practices that include the application of natural 
fertilizers, such as fertilizer made from fish parts, and 
exclude chemical fertilizers, herbicides, and insecticides. 



Before seedling transplanting in the middle of May, the 
natural fertilizers are applied to WF fields. Some of the 
fields that did not receive herbicide application were 
dominated by weeds, such as Monochoria vaginalis 
var. plantaginea. So farmers do the mechanical weeding 
once or twice in June. WF fields are flooded all 
year. Even in winter (not rice-growing season) WF 
fields are flooded for providing suitable habitats for 
waterfowl (Kurechi 2007). A typical farming schedule 
in the conventional paddies is as follows. Fungicides 
(Tetrachloroisophthalonitrile or hydroxyisoxazole and 
Metalaxyl) are applied to young rice plants in seedling 
boxes, and herbicides (indanofan, bensulfiironmethyl, and 
clomeprop (5 kg ha"^) or Pyriminobac-methyl, bromobu- 
tide, bensulfuronmethyl, and pentoxazone (10 kg ha '^)) are 
applied directiy to the field after seedling transplanting in 
early May. Insecticides (etophenprox (2-3 1 ha"^)) for rice 
leaf beetie Oulema oryzae Kuwayama are applied to the 
fields and then water is drained from the field for about 
ten days to stabilize rice yield and increase soil hard- 
ness in early June. After heading of rice plants in early 
August, insecticides for pecky rice bugs (dinotefuran 
(30 kg ha"^)) are applied once or twice to the 
field. The insecticides are applied to many of the con- 
ventional rice fields by radio-controlled helicopters. 
Most conventional paddies are seasonally dried-up 
after harvesting. 

Field survey 

A total of 31 WF and 7 conventional paddies were se- 
lected in an area of approximately 20 km^. The mean 
(SD) area of the study fields was 2100.8 ± 438.2 m^ and 
1780.3 ± 523.7 m^ for the WF and conventional fields, 
respectively. Rice seedlings were transplanted in May 
2007 at a mean (SD) density of 21.20 ± 2.0 hills/m^ and 
23.24 ± 2.0 hills/m^ in the WF and conventional fields, 
respectively. In August 2007, net samplings (20 sweeps 
with a 36-cm-diameter insect net) for spiders and their 
potential prey were performed at the center of each field 
just before and after using insecticides for pecky rice 
bugs in conventional fields (insecticide application: 6, 
18-19 August, first survey: 1-4 August; Second survey: 
15-19 August). The sweeping sampling is an excellent 
method to collect invertebrates inhabiting the upper- 
most vegetation layer in paddies easily. Invertebrate 
body mass was calculated as mass (mg) = 0.0305 x (body 
length (mm))^'^^ (Rogers et al. 1976) to measure the 
abundance of potential prey for spiders. As indices of 
habitat complexity, we measured rice height of six ran- 
domly chosen hOls and counted the number of individuals 
of the most abundant weed (M. vaginalis) within a 1 x 3 m^ 
quadrant at the center of each field (e.g., Rypstra et al. 
1999, Miyashita and Takada 2007). 
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Table 1 Independent variables used in the generalized 
linear mixed effect models for identifying limiting factors 
of spider density and species richness in winter-flooded 
(WF) paddies 



Independent variables 


Mean 


Range 


Tolerance values 


Habitat complexity 








Mean rice heiglnt (cm) 


67.34 


Dy.6/-//.83 


0.903 


Monochoria vaginalis 
(no.individuals/3 m^) 


67.42 


0-393 


0.794 


Prey availability 
(mg/20 sweeps with 
an insect net) 








CInironomid 


1 30.68 


14.23-975.86 


0.837 


Otiier nematocera 


32.07 


6.58-100.54 


0.913 


S. rubrovittatus 


7.08 


0-65.37 


0.746 


Brachycera 


5.23 


0.38-29.90 


0.831 


Hopper 


3.98 


0.79-10.17 


0.823 



paddy. The significance of each fixed effect was checked by 
using analysis of deviance and by checldng for the reduc- 
tion in residual deviance by reducing each effect by means 
of the Chi-squared test. 

The factors limiting the density and diversity of spiders 
in the WF paddies were analyzed by using generalized 
linear mixed effect models with Poisson error distribu- 
tion because the dependent variables are count data. We 
examined the limiting factors of spider density and spe- 
cies richness using only the data of the first survey to 
avoid redundant analyses, because spider density in the 
WF paddies tended to be higher in the first field survey 
than in the second (see Results) and the results obtained 
from preliminary analyses using the data of the second 
survey were similar to those of the first survey. This 
model included the number of individual spiders or 
spider species richness as the dependent variable, and 



Statistical analysis 

To evaluate the effects of the farming practices (WF vs. 
conventional) and the insecticides applied to conventional 
fields (before vs. after) on the numbers of individual spiders 
and spider species in paddies of both types, we applied gen- 
eralized linear mixed effect models with Poisson error dis- 
tribution because the dependent variables are count data. 
Fixed factors in the models were the main effects of the 
farming practice and the insecticide applications and their 
interaction, and the random factor was the effect of the 



Table 2 Mean density of each spider species or family per 
field for WF and conventional paddies 



Species/family 



Before insecticide 
application 



After insecticide 
application 





WF 


Conventional 


WF 


Conventional 


Tetragnatiia caudicuia 


2.774 


0.143 


0.406 


0 


T. praedonia 


0.774 


0.143 


0.469 


0 


T. extensa 


0.194 


0.143 


0.0406 


0 


T. maxiiiosa 


0.484 


0 


0.094 


0 


T. squamata 


0.097 


0 


0.031 


0 


Pacinygnatiia 
quadrimacuiata 


0.516 


0 


0.406 


0 


P. cierci<i 


0.161 


0 


0.344 


0 


Lycosidae 


0.065 


0 


0.250 


0 


Neoscona adianta 


0.065 


0 


0 


0 


Argiope bruennictiii 


0.032 


0.286 


0 


0 


Larinioides cornutus 


0 


0 


0.063 


0 


Linyphiidae 


0.032 


0 


0.281 


0 


Thomiside 


0.032 


0.143 


0.031 


0 


Theridiidae 


0 


0 


0.031 


0 


Unknown spiders 




0 


1.781 


0.143 
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Figure 1 Spider diversity in WF and conventional paddy fields. 

Boxplots for tlie spider density (a) and species ricliness (b) in the 
WF and conventional paddy fields before (filled bars) and after 
(open bars) applying insecticide to conventional paddies in 
early August. 
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two indices of habitat complexity (mean rice height and 
M. vaginalis density) and the biomass of five prey groups 
as fixed factors (Table 1), with the paddy field as a random 
effect. All tolerance values between the fixed factors were 
greater than the critical value of 0.1 (Table 1), indicating no 
significant coUinearity between them (Quinn and Keough 
2002). Model selection for the generalized linear mixed 
effect models was performed using Akaike's Information 
Criterion (AIC). AIC was calculated for each possible com- 
bination of fnced factors, and the models with small AIC 
were chosen as optimal models. All statistical analysis was 
carried out using R for Windows 2.13.1 (R Development 
Core Team 2011). 

Results 

A total of 234 spiders in 15 species or families were 
observed. Dominant spider species were Tetragnatha 
caudicula, T. praedonia, and T. extensa for web-building 
spiders and Pachygnatha quadrimaculata, P. clercki, and 
Lycosidae spp. for cursorial spiders (Table 2). Juveniles 
were identified only to family. 

The main effect of farming practice (WF vs. conven- 
tional) on spider density (x? = 23.369, P< 0.001) and 
species richness (xi = 17.817, P < 0.001) was highly sig- 
nificant, while the interaction effect of farming practice 
and insecticide application was significant on neither 
spider density nor species richness (density: Xi = 2.792, 
P = 0.095; species richness: Xi = 3.204, P = 0.073). This 
result indicates that spider density and species richness 
were higher in WF paddies than in conventional 
ones both before and after insecticide applications to 
conventional paddies (Figure 1). Since the main effect of 



insecticide application (before vs. after) was significant 
only for spider density (density: Xi = 5.920, P = 0.015; 
species richness: Xi = 1-417, P = 0.234), insecticide appli- 
cation to conventional paddies was suggested to reduce 
spider density in both the WF and conventional paddies 
(Figure 1). Moreover, the variations in spider density and 
species richness in the WF paddies seemed to increase 
after applying insecticide: the interquartile range — the 
difference between upper (75% quartile) and lower (25% 
quartile) hinges — of spider density was 5.5 before in- 
secticide application, but 6 after (Figure la), and the 
interquartile range of spider species richness was 2 be- 
fore insecticide application, but 4 after insecticide appli- 
cation (Figure lb). 

The best model for explaining the spider density and spe- 
cies richness in WF paddies included the biomass of chi- 
ronomids and other nematocera (Table 3) but did not 
include any factors of habitat complexity, such as mean rice 
height or M. vaginalis density. These two groups of prey in- 
creased both spider density and species richness (Figure 2). 



Discussion 

We revealed that spider density and species richness 
were higher in the WF paddies than in the conventional 
ones both before and after insecticide application to 
conventional paddies for pecky rice bugs in early August. 
Moreover, our results suggest that the insecticide appli- 
cation decreased spider density in not only conventional 
paddies but also WF ones. The spider diversity in the 
WF paddies was significantly affected by the availability 
of two prey groups, chironomids and other nematocera. 



Table 3 Estimated coefficients(SE) and information theoretic statistics of the top five and null (intercept only) models 
explaining spider density and species richness in WF paddies 



Rank Intercept 


Rice 
height 


M. vaginalis 


Chironomid 


Other S. rubrovittatus 
nematocera 


Brachycera 


Hopper 


AIC 


A AIC 




Spider density 




















1 -0.941 (0.457) 






0.473 (0.092) 


0.014 (0.004) 






52.10 


0 


0.119 


2 -0.717(0.471) 






0.460 (0.090) 


0.013(0.004) 


-0.028 (0.022) 




52.45 


0.34 


0.100 


3 0.844 (1 .736) - 


■0.027 (0.26) 




0.487 (0.93) 


0.013(0.004) 






52.97 


0.87 


0.077 


4 -0.755 (0.543) 




-0.031 (0.53) 


0.445 (0.0102) 


0.014(0.004) 






53.75 


1.65 


0.052 


5 






0.473 (0.091) 


0.013(0.004) 


-0.025 (0.023) 




53.80 


1.70 


0.051 


Null 1.490(0.164) 














73.44 


21.33 


0 


Spider species richness 




















1 -0.535 (0.481) 






0.278 (0.098) 


0.008 (0.005) 






21.13 


0 


0.074 


2 -0.286 (0.517) 






0.262 (0.097) 


0.009 (0.005) 




-0.053 (0.044) 


21.63 


0.50 


0.057 


3 -0325 (0.451) 






0.294 (0.096) 








21.87 


0.74 


0.051 


4 -0372 (0.515) 






0.268 (0.099) 


0.007 (0.005) 


-0.020 (0.25) 




2240 


1.27 


0.039 


5 -0.138(0.481) 






0.278 (0.097) 




-0.025 (0.025) 




22.73 


1.60 


0.033 


Null 0.960(0.111) 














28.68 


7.55 


0.002 



Ranking of the sub models is based on Al<aike's information criterion. 
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Figure 2 Spider diversity and their potential prey in WF paddy 
fields. Relationships between total biomass of important prey 
(chironomids and otiier nematocera) and spider density (a) and 
species riciiness (b) in WF paddies before applying insecticide to 
conventional fields. Poisson regression lines are also shown. 



The significant main effect of farming practice (WF vs. 
conventional) on spider density and species richness 
appears to indicate that these values were already lower 
in conventional fields than in WF ones before the in- 
secticide application to conventional paddies (Table 2, 
Figure 1). At least four possible explanations for the 
farming practice effect can be identified. First, fungicides 
applied to young rice plants in seedling boxes, herbicides 
and insecticides early in the rice-growing season may 
be responsible for the reduced spider diversity in con- 
ventional fields during the early-growing season. Such 
agricultural chemicals are known to affect invertebrates 
(Oyama and Kidokoro 2003; Wilson et al. 2008; Amano 
et al. 2011) and were probably responsible for decreas- 
ing the number of spider prey in the paddies. The 
second mechanism for the farming practice effect on the 
spiders is due to the drainage of conventional paddies in 



early June. Drainage dries and hardens the soil in con- 
ventional paddies (see Study site) and likely decreases 
the number of aquatic invertebrate, including juveniles of 
dipterous insects, which are important prey for spiders 
(Ishijima et al. 2006; Tahir and Butt 2009). Third, organic 
fertilizer applied to the WF fields before seedling trans- 
planting might increase spider diversity due to increasing 
the availability of important prey for spiders. Previous stud- 
ies have reported that detrital subsidy, such as organic 
fertilizer, enhances spider populations indirectly through in- 
creasing detritivores, such as chironomids and coUembolas 
(Murata 1995; Setde et al. 1996; Halaj and Wise 2002). 
Finally, tubificid worms as well as waterfowl fertilize WF 
paddy soil with their feces (Kurechi 2007), which may in- 
crease spider diversity due to increasing the availability of 
important prey for spiders. Because WF fields are flooded 
all year, tubificid worms are known to active even in 
winter, fertilizing the soil in WF fields. On the other 
hand, conventional fields are dried-up after harvesting. 
Clarification of these causalities requires manipulative 
experiments controlling each component with many 
replications. 

The main effect of insecticide application (before vs. 
after) on spider density was also significant, suggesting 
the insecticides applied to conventional fields may de- 
crease spider density in not only conventional but also 
WF paddies, despite not applying insecticides to WF 
paddies directly. One possible explanation is as follows. 
WF and conventional paddy fields are spatially inter- 
mingled in the study area, so insecticides for pecky rice 
bugs sprayed onto the conventional rice fields by radio- 
controlled helicopters may easily diffuse into adjacent WF 
paddies and consequently affect spider density. Increased 
variation in the spider density and species richness in WF 
paddies after insecticide application (Figure 1) appears to 
support this possibility, because the intensity of the negative 
effect of the insecticides on spiders inhabiting the WF pad- 
dies may vary depending on the distance between the WF 
and conventional paddies. If so, increasing the spatial 
gathering of WF paddies will lead to an increase in the 
number of spiders in each WF paddy due to avoiding 
the diffusion of insecticides sprayed by radio-controlled 
helicopters. Enhancement of spider populations by stra- 
tegic spatial placement of WF paddies may also suppress 
damage by pecky rice bugs in WF paddies. Spider spe- 
cies with high abundance in the WF paddies were also 
shown to predate pecky rice bugs through previous 
DNA-based gut-content analysis (Kobayashi et al. 
2011; Takada et al. 2013). Moreover, Tetragnatha spp. 
spiders were revealed to suppress bug density and the 
consequent damage in organic paddies fields (Takada 
et al. 2012). 

Spider diversity in WF paddies increased with 
the increasing abundance of chironomids and other 
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nematocera. This was revealed by statistical analysis of 
the limiting factors of spider diversity (Table 3). Be- 
cause we used a comparative approach, the casual rela- 
tionship between spiders and these insects was no 
proven in this study. However, previous studies have 
suggested that these dipterous insects are frequently 
consumed by the spiders in paddies (Ishijima et al. 
2006; Tahir and Butt 2009). Taken together, the dipter- 
ous insects most likely limit spider diversity in WF 
paddy fields. Our study is one of the few that provides 
evidence of the correlation between the abundance of 
important prey and increased spider diversity in pad- 
dies, with many replications. 

Spiders and their potential prey were collected only by 
one method (i.e. sweeping sampling) in this study. 
Sweeping sampling is an excellent method to collect in- 
vertebrates inhabiting the uppermost vegetation layer in 
paddies easily. But it is necessary to use other methods 
(e.g., direct counting), to estimate spiders inhabiting the 
bottom layer, such as Lycosidae. 

Many studies have shown that flooding rice paddies 
can provide suitable habitats for waterfowl (Elphick 
and Oring 1998; Tourenq et al. 2003; Kurechi 2007). 
The present study suggests that in WF paddies the di- 
versity of generalist predators is higher than in con- 
ventional paddies during the rice-growing season. 
Moreover, a combination of management at both the 
landscape and field level seems to be more effective 
for increasing spider diversity in WF paddies. Spiders 
are not only important natural enemies of various in- 
sect pests (Kiritani et al. 1970; Oraze et al. 1989; 
Takada et al. 2012), but also important prey for larger 
predators, such as frogs (e.g., Hirai and Matsui 2001) 
and larger spiders (e.g., Ishijima et al. 2006) in paddies. 
Taken together, the enhancement of spider diversity 
in WF paddy fields might achieve both biological con- 
trol and biodiversity conservation in rice paddy 
ecosystems. 
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